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FIELD OF THE INVENTION 

The present Invention relates to a method for channel estimation according to the 
preamble of claim 1. 

Such a method may be used In particular In any communication system using UMTS 
standard* 

BACKGROUND OF THE INVENTION 

A UMTS communication system, defined in the UMTS standard, comprises at least a 
base station comprising a transmitter that sends a signal to users, that is to say to receivers 
through downlink communication channels, and receivers that send a signal to said 
transmitter of said base station through uplink communication channels, 

A downlink communication channel uses two logical channels, a first channel called 
common pilot channel (CPICH) that continuously transmits known pilot symbols in particular 
for controlling the transmission of a signal and a second channel called dedicated physical 
channel (DPCH). Said second channel DPCH comprises two dedicated channels, a first sub- 
channel called dedicated physical data channel (DPDCH) time multiplex with a second sub- 
channel called dedicated physical control channel (DPCCH). 

Transmit beamformlng can be employed at the base station to Improve the second 
logical channel DPCH transmission. It consists in providing dedicated beams for useful data 
users. It Improves the capacity of the UMTS communication system to manage a certain 
number of users by reducing interferences between twd users: On the contrary, 
beamformlng is not suited for the first logical channel CPICH transmission that is to be 
uniformly broadcast In time and space to all users. As seen from the receiver, the second 
logical channel DPCH is in general different from the first one CPICH. Usually, when no 
transmit beamformlng is used, in order to optimize the receiver, one uses a maximum a 
posteriori probability (MAP) optimization, i.e. one uses channel estimates of the CPICH 
channel as channel estimates for the DPDCH channel. But, when transmit beamformlng is 
applied, channel estimates obtained on the CPICH cannot be directly applied as channel . 
estimates for the DPCH as it is usually done because they both differ. 

SUMMARY OF THE INVENTION 

Accordingly, it is an object of the invention to provide a method for channel 
estimation with transmit beamforming, said channel comprising a first and second channels, 
which achieves an efficient optimization of a receiver without too much complexity. 

To this end, there is provided a" method as claimed in claim 1. 



As we will see in detail further on, such a method enables to optimize a receiver, i.e. 
reducing interferences between users, by calculating an estimate of the first physical 
channel from the second physical estimate and moreover from a first logical channel 
estimate. The use of a first logical channel estimate allows a better tracking of the channels 
variation In time, said first logical channel transmitting continuously contrary to the second 
physical channel. 



BRIEF DESCRIPTION OF THE DRAWINGS 

The invention and additional features, which may be optionally used to Implement 
the invention to advantage, are apparent from and elucidated with reference to the 
drawings described hereinafter. 

- Fig. 1 represents a diagram illustrating the main steps of the method for channel 
estimation according to the invention, and 

- Fig. 2 illustrates an estimator using the method for channel estimation according to the 
invention. 

DETAILED DESCRIPTION OF THE EMBODIMENTS 

In the following description, well-known functions or constructions by a person of 
ordinary skilled in the art are not described in detail since they would obscure the invention 
in unnecessary detail. 

The present invention relates to a method for channel estimation with transmit 
beamforming, said channel comprising a first CPICH and second DPCH channels. Said 
method is used In particular In a RAKE receiver of a UMTS communication system, for 
example in a receiver of a mobile terminal MT. 

A UMTS communication system comprises at least a base station BS comprising a 
transmitter that transmits a signal to users, that is to say to some receivers through 
downlink communication channels, said receivers sending a signal to said transmitter of said 
base station BS through uplink communication channels. The transmitted signal comprises 
useful data and pilot symbols. 

A downlink communication channel uses two logical channels, the first one is called 
common pilot control channel CPICH that continuously transmits known pilot symbols and 
that employs a fixed spreading factor M equal to 256, and the second one is called dedicated 



physical channel DPCH. Said second channel DPCH comprises with useful data user forming 
a first sub-channel called dedicated physical data channel DPDCH which are time multiplex 
within a same time slot with known pilot symbols forming a second sub-channel called 
dedicated physical control channel DPCCH. 

Transmit beamforming Is employed at the base station BS in order to Improve the 
DPCH reception at the RAKE receiver of a mobile terminal MT. It consists in providing 
dedicated beams for useful data. It improves the capacity of the UMTS communication 
system to manage a certain number of users by reducing interferences between two users, 

A RAKE receiver comprises fingers that are associated with each path of a beam 
dedicated to useful data that has enough power, a finger being a filter adapted to recover a 
path. 

It will be noted that logical channels use physical channels in order to transmit their 
data. Classically a physical channel has a temporal structure and a spatial structure, a 
physical channel being for example an antenna. 

It will be also noted that the transmitted signal is composed with bits, said bits being 
first encoded, the encoded bits being called symbols. Then, these symbols are modulated 
(e.g. by a QPSK modulation, well know by the person skilled in the art). Then, said 
modulation symbols, are submitted to a spreading. The use of a spreading factor M and of 
an associated code C enable to transmit a signal on a larger bandwidth. After the spreading, 
the bits are defined in term of a chip. We have the following formula: bit rate * M = chip 
rate. 

In order to recover the useful data that have been received by the UMTS RAKE 
receiver, one has to estimate the transmission channels such as the downlink user dedicated 
physical channel DPCH. While such an estimate is done, effects on the received signal due to 
Interferences of such channels can be removed in order to recover said useful data. 

As we will see, for such an estimate, there is a first step of model of the physical 
channels associated said first CPICH and second DPCH channels, raid models being a linear 
superposition of discrete multipath MP components (characterized by the term 
Cp(t)$(r-T p ), as we will see below) following an uncorrelated-scattering wide-sense 

stationary model. A multipath component is characterized by a time-varying complex 
coefficient c p (t) and a propagation delay r p . Said time-varying complex coefficient has a 

CPICH multipath coefficient c ich (t) and a p-th DPCH multipath coefficient 



Cpdpchi*) = P ? c pcpi<*Sf) wh,ch depends on a complex factor fi p due to the transmit 
beamforming of the base station BS. 

Signal Model 

We first begin to model a general channel h. This general model will be applied to 
the CPICH and DPCH channels afterwards. 

1) General Sjgnql Model 

We consider a general case of the transmission of a signal 5(0 through a muitJpath 
channel h. A continuous-time complex base-band received signal at a RAKE receiver is given 
by 

y(0 = *(0+v(0 & 



v(0 defines noise plus interference, 

y(t) Is a received signal before despreadlng and before an A/D conversion, 
x(t) is a received signal that comprises the useful data, which have been subjected 
to a distortion due to the effects of the channel h. 



is given by 

*(0= J*(M-T)z(r)dr ( 2 ) 

where /<£t) is a time-varying mulUpath channel impulse response and 

*(0 = £«(*M<-*r c ) O) 



where: 

- z(t) is a signal just after a D/A conversion in the transmitter, 

- T e , is a chip period, 

- a(k) = s([k/M J)* d(k), M denotes the spreading factor, and the term 

L J is a rounded down, 

- 5(/7) is an n-th modulation symbol (e.g., QPSK), 



- o(A) Is a *-th chip, and 

- ytQ is a pulse-shaping filter that limits the base-band of the receiver. 
Note that this model can be used by systems with several spreading layers, like DS- 

CDMA system defined by the standard IS-95. 

The channel h is assumed to follow a wide-sense stationary uncorrected scattering 
model with Rayleigh fading and multlpath response as described in the document hereby 
incorporated by referenced. G. Proakis Digital Communications, NY: McGraw-Hill, 3rd 
edition, 1995" which is incorporated by reference. It means that: 

Kt 9 T) = %C p (t)S{T-T p ) (4) 
/>=! 

where P denotes a number of paths available for one user, and qHQ represents a 

path time-varying complex coefficient at propagation delay tp> The path complex 

coefficients c^O vary depending on the velocity of the mobile terminal MT and on 

p 

surrounding scattering objects such as buildings, mountains, cars, trains.... 2 IqXOI 2 

corresponds to the energy that has gone through a DPCH channel. These complex 
coefficients are not known by the receiver, because the physical channel characteristics 
such as the complex coefficients are unknown from said receiver. Therefore, they have to be 
estimated as we will see later on. 

In the case of transmit beamforming through a Q-antennas array of a RAKE 
receiver, Q being the number of antennas of the array also called rake, the channel model 
expressed by (4) becomes 

Kt,r) = £w H a(9 p )c p (t)S{T-T^ (5) 
where: 

- h> = O 0 , . . . , w Q _ x f denotes a transmit beamforming weight vector, which enables to 
weight the antennas when forming the beams, 

- a(0)=[a o (0),... > a fi -i(0)] r denotes an antenna array response vector in a direction 6, 

and 



- The superscripts ( . ) r and ( . )" denotes transpose and Hermitian transpose 
respectively. 

By Inspecting equations (4) and (5) one may notice that, during a time slot, transmit 
beamforming only affects the amplitude of the path time-varying complex coefficients CfH§ 

5 by a complex factor w* a(0 p ) and not the delays It means that the channel's temporal 

structure is not affected by the antenna, neither by the transmit beamforming. 

Moreover the path time-varying complex coefficient q£f) stays the same in both 
case (3) and (4) up to a complex factor^ = w H a(0 p ) . 

10 2) Received Siqp al Model for DPCH and CPICH 

We now apply this general channel model to the DPCH and CPICH channels as 

follows. 

We consider the case where the first CPICH and second DPCH channels are 
15 superposed, that is to say only when the useful data are transmitted, knowing that the first 
CPICH channel continuously transmits pilot symbols while the second DPCH channel only 
transmits when useful data have to be transmitted. 

In that case, the continuous-time complex base-band received signal (1) can be 
written as 

20 y(t) = Xcpieh (0 + x dpch (0 + v(/) (6) 

where Ajp/wCfl and Ai^(0 are the signal components associated with the first 
CPICH and second DPCH channels respectively given by 

25 where: 
and 

k 

k 



where: 



- W*)"WL*'^«J>^W*)' 

- Scpain), Sjpctiri) are the /rth modulation symbols, assodated with the first CPICH and 
second DPCH channels respectively, and 

- dcptdJil , ddpdify are the *-th chips, associated with the first CPICH and second DPCH 
channels respectively. 

It will be noted that the spreading factor M may not be the same for the first CPICH 
and second DPCH channels. In the UMTS standard the first channel CPICH has a fixed 
spreading factor equal to 256 whereas the second channel DPCH has a spreading factor that 
can vary between all powers of 2 between 4 and 512. 

Furthermore, one can notice that the CPICH is not modulated, i.e., s^(/7)=yl(l+j) 
for any n symbol, where A Is an amplitude normalization factor. It means that s^ri) 
doesntvary in time. 

We assume now that the RAKE receiver has perfectly recovered the path delays as 
the RAKE receiver knows the delays r p of the paths associated to the CPICH and DPCH 

channels. The RAKE receiver comprises an analog/digital converter that allows to go from a 
continuous time (t) to a discrete time (n). 

We also assume that the RAKE receiver has already de-spread the received signal 
y(t) in order to detect the CPICH and the sub-channels DPCCH/DPDCH channels. It means 
that the RAKE receiver has multiplied the received signal y(tpy the code C associated with 
the spreading factor M used for the transmission of a signal so as to cancel the effect of the 
spreading and so as to restore the original narrow bandwidth. During the dispreading, a 
sampling Is done according to a symbol rate 1/T, T representing the symbol period. 

For the sake of simplicity, we limit the treatment to a single path/finger case of the 
receiver. Thus, we will drop the subscript p to denote the time-varying complex coefficient 
hereinafter also called path fading channel coefficient. 

Furthermore, in a first Instance, we limit the analysis to the case where the second 
channel's DPCH spreading factor Mis equal to the one of the CPICH. We shall generalize the 
analysis to other spreading factors successively. 

Assuming to model all interferences and especially multiple-access interferences 
(MAI) and inter-path Interferences as If it was additive white Gaussian noise (AW6N), the 
de-spread signals, also called observations for a RAKE finger related to the CPICH and the 
DPCH channels are respectively given by 
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y cpl cH(n) = c(n)s cplch (_n) + v eptch (n) 
^W = w"a(e)c(n)5 4lcl (n)+v iM (») 
where: 

- v^/?) and veptetin) denote noise-plus interference terms associated with the second 
DPCH and first CPICH channels respectively, classically assumed to be independent 
from each other, 

- Is the path time-varying complex coefficient, 

- Sapain) and denote the pilot symbols of the DPCH and CPICH channels 
respectively, at a symbol time n. 

It will be noted that a MAI is for example interferences between different users or 
interferences for one user due to the first channel CPICH or due to different data flows 
dedicated to this user, and that an inter-path interference is due to interferences between 
path impulse responses at different delays r . 

Under the assumption that the first CPICH and second DPCH channels vary every 
symbol time n and that the interference and the noise can be modeled as an additive white 
Gaussian noise, we then make a maximum likelihood (ML) static channel estimates at every 
symbol time n. This ML can be built as follows: 

Ccptch (») = <W 00>VA (») = C <») + V cpich (») (8) 

00 = s lc (n)y dpch («) = yw " a(fl)c(n) + v dpch («) 

where, v^OO and v^Cn) symbolize the AWGN associated with the first CPICH 

and second DPCH channels respectively, and where without loss of generality, we assumed 
the pilot symbols of both channels to be normalized such that |s$k<</?)I = \Scptdiri)\ = 1. 

Note that s*(n) is a pilot symbol multiplied by Its complex conjugated. 

Note that the pilot symbols a**(?i) and SdpchQt) of both channels CPICH and DPCH 
are known by the transmitter and RAKE the receiver, y is a coefficient that symbolizes a 
difference of energy between the pilot symbols of the second sub-channel DPCCH channel 
and the useful data of the first sub-channel DPDCH channel. 

Upon the models of the first DPCH and second CPICH channels, the following 
remarks can be made. 

1. Such ML static channel estimates are formed by demodulating 



(multiplication by the complex conjugated) observations and at the output 
of a finger in data aided (DA). But it can be done for decision direct (DD) fashion. 
While DA demodulation is achieved using the known pilot symbols provided on the 
second sub-channel DPCCH and first logical CPICH channels, DD demodulation 
requires the use of hard decisions taken at the RAKE receiver output as pilot 
symbols, when known pilot symbols are not available, as it occurs during the first 
sub-channel DPDCH channel transmission. It means that when receiving some bits 
of a signal, a hard decision is taken in order to define a symbol. Such symbol is 
taken as a pilot symbol. 
Therefore, the formula (8) can be written as: 

s dpck (j*)y<ipch ( w ) ~ ^(n)y dpcch (n) +S*^(n)y dpdch (ri) 
= wV0)c(w) + v^ A (n)+^ 

2. During a symbol rate 1/T, the path time-varying complex coefficients of the first 
CPICH and second DPCH channels only differ by a complex scalar factor given 
by/^yw^fl). 

3. When there Is a transmission, there is a transmit power offset between both first 
CPICH and second DPCH channels because the power used to transmit a signal is different 
for said first CPICH and said second DPCH channels. While the power is relatively constant 
for the first channel CPICH, the power for the second channel DPCH is changing as a 
function of some propagation conditions. For example, when the mobile terminal MT is near 
to the base station BS, the power is lower than when said mobile terminal MT is far from 
said base station BS. As said power offset between CPICH and DPCH is a real and is 
constant over a time slot, It can be taken into account in the complex factor p. Hence, DPCH 
power assignment is included by properly normalizing the term 

4. On a physical point of view, to each path on the first channel CPICH corresponds 
a path on the second channel DPCH with the same delay r and a same Doppler spread but 
with different complex paths' amplitudes p and p c(f) and so p c(n), the Doppler spread 
being related to the Doppler effect made by the mobile terminal's MT movement or made by 
the scattering objects in movement around Hie mobile terminal MT. The two complex paths' 
amplitudes are only different up to a complex factor jffthat is determined by the scalar 
product of the beamforming weight vector ivand the antenna array response a(0) in 
correspondence of the direction 6 of that path. It means that the two channels CPICH and 
DPCH behave In the same way up to said complex factor p 9 which is modified as a function 
of a path. Hence, these channels CPICH and DPCH are not different in a temporal point of 
view. 
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5. On a mathematical point of view, it means that the path time-varying complex 
coefficients and so c(n) vary as a function of the time as it has been mentioned before, 
as said coefficient depends on said Doppler spread. 

5 From these remarks, we can conclude that, as the two channels CPICH and DPCH 

behave in the same way up to said complex factor p, said complex factor p can be used to 
estimate more particularly the first physical channel PPDCH. Hence, we have to calculate 
first an estimate of said complex factor p. 

10 To tills end, we use the DPCH pilot symbols jointly with the pilot symbols provided 

by the CPICH for estimating the unknown complex factor p. Once an estimate p of said 
complex factor p is calculated, a channel estimate for the first sub-channel DPDCH channel 
will be easily obtained by multiplying the estimate p of said complex factor by the channel 
estimates obtained over the CPICH channel, 

15 

It will be noted that this approach can include the DA ML only or with the DD ML 
channel estimates according to equation (8) in the estimate of p. 

Note that this approach is more advantageous than the one that consists in using 
20 only the pilot symbols provided by the second sub-channel DPCCH to directly estimate the 
first sub-channel DPDCH. Indeed, the number of pilot symbols provided by said second sub- 
channel DPCCH is usually fairly small (from 1 to a maximum of 16 depending on the slot 
format), and in some cases not enough to provide sufficient noise suppression for a reliable 
DPDCH first sub-channel estimate. Moreover, as the second sub-channel DPCCH is not 
25 continuously transmitting pilot symbols, at high mobile terminal T velocity, there is a risk of 
not taking Into account some variations of the first sub-channel DPDCH during its 
transmission while said second sub-channel DPCCH is not transmitted, such that the channel 
estimates corresponding to said sub-channel bPCCH transmitting period cannot be applied 
as channel estimates for the whole successive first sub-channel DPDCH period transmission. 
30 With the first channel CPICH that transmits continuously some pilot symbols, those 
drawbacks are overturned. 

DPDCH Channel estimate from CPICH and DPCCH Channel estimates 

35 Therefore in order to estimate the first sub-channel DPDCH, we process each 

multipath component as follows. 
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1) In a first step, we build DPCH and CPICH instantaneous maximum likelihood ML 
channel multipath MP coefficients estimates c dpch {n), c cpich (n). 

As above described in the paragraph Received Signal Model for DPCH and CPICH, 
the estimates are: 

Ccpick in) = {n) ycpich (n) = c(n) + v cpich (n) 

3 

Cdpch (») = J^rf (n)y dpeh (w) = yw* a(0)c(n) + v dpeh («) 

In a preferred embodiment, we use the DA fashion: The DPCH Instantaneous 

maximum likelihood ML channel multipath MP coefficient estimate c dpch (n) is only based on 

the ML coefficient estimate of the second sub-channel DPCCH. That is to say, only the 
estimate of the second sub-channel DPCCH is used jointly with the estimate of the CPICH 
10 channel. 

Hence, the estimates are: 

Ccpich in) = s] pich OOy^ in) = c{n) +y cpich (n) 

Zdpch in) = in)y dpch (n) = yv H a(0)c(n) + v dpch (n) 

2) In a second step, we build an estimate c^ JcA (n) of the CPICH multipath 

15 coefficient 

In a non-limitative preferred embodiment, the CPICH multipath coefficients estimate 
is done by a linear prediction as described in the document hereby incorporated by 
reference [1] W J. Bsltersee, G. Fock, P. Schultz-Rittich, and H. Meyr, "Performance analysis of 
phasor estimation algorithms for FDD-UMTS RAKE receiver," IEEE (P Symp. on Spread 
20 Spectrum Teehonofogfes and Applications, NJIT, NJ, September 2000." 

Under the assumption that the channel paths are Gaussian complex random 
process, a maximum a posteriori probability (MAP) channel path coefficient estimate is given 
by a conditional expected c mp (n) value given the whole observation y(n) 

Yin) = [ycpich i n ~ K c P ich + 0 - ycpichin) ] r and the corresponding pilot 

25 symbols sequence $(n) = [s^ (n - K cp(ch +1) ... (n) ] T such as 

c mp (n) = E{c ML (n) | y(n), s(n)} (9) 

where c^in) represents the Instantaneous ML static channel estimate, is the 
30 number of available pilot symbols from the CPICH channel. 
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In the document referenced [1], It Is shown that a MAP optimal channel estimate (9) 
can be built by low-pass filtering instantaneous ML channel estimates (8), i.e., 
Cmp^f"c ML {n) (10) 

wherec J1/i («) = [c (i)fcft (») ... (n - + 1) fdenotesthe 

vector of the ML channel estimates and f=[f(N-V) ... /(0)f, with N^K^, denotes 
the low pass filter coefficients respectively. 

Although the current formulation corresponds to a classical forward linear prediction 
(we go from n to n-K+1...), we shall remark that it can trivially extent to both, backward or 
forward-backward (i.e., two-sided) linear prediction. 

The linear prediction is well known by the person of ordinary skill in the art and 
therefore wont be further described. 

Note that for a DD fashion, we will have the formula (9) extents to the DPCH 
channel such as: 

y(")=[y <pfc *(»-** te »+ 1 ) ••• ycptchW ^(»-^*c» +1 > - y^ n )f 

and the corresponding pilot symbols 
sequence 

s(") = [^ W») ^{n-K^+l) - s^in)} 7 such as 
Cmap(») = EtftaWyinXsin)) (9) 

where c ML {n) represents the Instantaneous ML static channel estimate, K^, lch and 
K dpch is the number of available pilot symbols from the CPICH channel and DPCCH sub- 
channel respectively. 

3) In a third step, we build estimates of a DPCH and CPICH cross-correlation 
E{c 4>eA (»)c qridi (n -/)*} and of a CPICH auto-correlation E^^OOc^Oi -/)*} from 
maximum likelihood ML estimates of step 2 at non-zero lag for noise suppression (t*0). 

The estimates are built as follows. 

From the formula (8), we have: 
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Ccplch («) = C(«) + («) 

<W») = r«»mc(n) + y dpd ,in) ^ C ° mpte ^ faCtor giVen 

In a preferred embodiment, only the DA fashion is used: the cross-correlation of the 
DPCH and CPICH channels is a cross-correlation of the second sub-channel DPCCH and first 
channel CPICH. 

Hence, we have: 

a fr,\ - ™"«/zjw ^ with the complex scalar factor given 

by/? = w"a(0). 

The cross-correlation faik) of the DPCCH and CPICH ML channel estimates, at lag /, 

Is: 

I 

= E {(/} c(n) + V4**(w))(c («-/)*+ vcpuwn -/))} = 

E{0 c(n)c(n-l)* +V4 PC/ ,(n)*c(n-iy +J3 c(n)*vcpuh(n-l)+vdpch(n)»vcpmn-l)} 
One can notice that when we take the AWGN model, noise samples v^ lch (») and 
v dpch (") taken at any lag l* o are uncorrected. Therefore, their means are null. Hence, 
= B{fic{n)c{n-l)*} 
= /3B{c(n)c(n-l)'} 
Thus, 0 = Ml)/B{c(n)c(n-r)*} 

The auto-correlation ^ (/) of the CPICH channel estimate, at lag /, is: 

#«(0- E{^(»)^ /cA (»-/)*} 

= E {(c(/i) + v«^(»))(c («-/)*+ vvmn - /))} = 

Since noise samples v^Qi) taken at any lag /*oare uncorrelated, their means 
are null. Hence, 

<U/) = E{c(»)c(«-/)*} 
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4) In a fourth step, we estimate said complex factor >0 with said estimates In step 3. 

We assume in a first instance that the DPCH and CPICH channels have the same 
symbol rate 1/T, i.e. the same spreading factor M. We shall extend the treatment to the 
more general case of different spreading factors successively. 

Consider the instantaneous ML channel estimates provided by (8). Then, under the 
assumption of modeling all the Interference and noise terms as an AWGN, the complex 
factor p can be estimated as follows. 



V/*0(11) 

^^("-0} feci 1 ) r ccO 

10 since again with the AWGN noise model, noise samples taken at any lag / * 0 are 

uncorreiated, the mean of such noise samples are equal to 0. Hence, one advantage is that 
we get rid of the noise and interference that annoy us on the CPICH and DPCH channels. 

In practice the true correlations r & (/) and r.(Z) cant be calculated, as a mean has 
15 to be calculated on an infinite number of samples. Therefore, these correlations need to be 
replaced by correlation estimates r^if) and r^if) . 

In a preferred non-limitative embodiment, we use sample moments as r ic (l) and 

r rf c(0 = ^-S^(")^(»-0 ' (12) 

20 and 



N p „«o 



r cc (/) = ^^.(^U(^ / ) (13) 



p 

so that 



N „ n=0 



where n p is the number of pilot symbols comprised in the second sub-channel 

25 DPCCH. 

Note that if the quality of the estimation is not good enough because of a too low 
value of the number of symbols at, , one may increase n„ by including DD channel 
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estimates, when reliable symbol decisions can be taken at the output of the RAKE receiver, 
that Is to say when there is a good signal to noise ratio SNR. 



In a preferred non-limitative embodiment, in order to improve the estimates, further 
noise rejection can be achieved by considering the weighed sum 

P-^LYiPi . L being an integer ♦ ( 15 j 

as an estimate for fi, where they, are weighed coefficients and are to be chosen In 
order to account for the reliability of the estimates p, . To derive the optimal weighed 
coefficients J*-,, one should determine the joint statistic of fi„ but that cannot be done in 
analytical way, as it would be too complex. 

But, one can notice that the quality of the estimates 0, is likely to decrease as the 
lag /Increases. Indeed, when the lag /increases, any channel correlation decreases whereas 
the noise Increases. Hence, for a fixed number of pilot symbols*, the noise bias in the 
sample moments estimates (9) and (10) increases as a channel correlation decreases. 

Thus, in order to limit the algorithm complexity, In a preferred embodiment, we can 
limit the estimator (15) to the first term, that is to say we regard the first term as having the 
optimal weighed coefficient y x = 1. Indeed, It would be too complicated to make a 
weighted mean on the other terms /},, 1*1. Hence, we have: 



P ~ 7F3 (16) 

££«u(")3U*(»-l) 



5) In a fifth step, from steps 2 and 4, we build a first sub-channel DPDCH mulUpath (MP) 
coefficient estimate c dpdch (n) = fi c^n) . 

In order to build the first sub-channel DPDCH estimate, for all the muWpath 
component MP, we repeat the steps 1 to 5. 
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poncH Chann e l Estimate at different symbol rates from CPICH 

and DPCCH 

As we have seen before, we have assumed in a first instance the first CPICH and 
second DPCH channels to have the same symbol rate 1/T. But, in reality, it is not always the 
case. When the second channel's DPCH spreading factor A/is not equal to 256, the symbol 
rates of the second channel DPCH and of the first channel CPICH differ. 

Hence, in a preferred non limitative-embodiment, the channel estimates of either the 
first channel CPICH or the second channel DPCH are interpolated to match the symbol rate 
of the other channel estimates in order to apply the previous algorithm (16) so as to 
estimate the complex factor p. 

This step of interpolation is done after the first step of the ML estimates of the DPCH 

and CPICH c^Cn), c^in) described before. 

In deriving the ML-statlc channel estimates (8), both first CPICH and second DPCH 
channels are implicitly assumed constant over the longer symbols period between said 
channels, i.e. between two cfnjand c(n+l). 

Therefore, in a preferred non-limitative embodiment, it is easier to make an 
interpolation to the lower symbol rate 1/T between said first CPICH and second DPCH 
channels. In particular, when the DPCH spreading factor A/is smaller than 256, the 
interpolation to the CPICH symbol rate Is obtained by averaging DPCH ML static channel 
estimates over a number of symbols equal to 256/A/. Whereas, when A*=512, the DPCH 
symbol rate being lower than the one of the CPICH, then the average Is taken on the CPICH 
over two symbol channel estimates. 

In another embodiment, interpolation can be done to the highest symbol rate 1/T 

between DPCH and CPICH. 

One can notice that these embodiments have the advantage of not coloring the 
estimation noise i.e. it doesnt imply change on the estimate algorithm (16), as we have 
seen that noise samples are uncorrelated. 

Estimator. Properties 

As we have seen, an estimator has been derived applying the well-known method of 
samples moments, where the cross- and autorcorrelation were replaced by the 
corresponding sample estimates. Moreover, the least complexity form of the estimator 
reduces to computing the cross- and the auto-correlation at one lag /as in the algorithm 
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(16). Unfortunately, even in this case statistical properties of the estimator cannot be 
derived in analytical way. It means that from the algorithm (16), we can not derived the 
performances of the estimator such as the number of pilot symbols Afcwe need to have to 
reach a precision of 1%, the average precision of said estimator.... However it is 
straightforward to see that the estimator is asymptotically unbiased. For AV*», 
4(0->'' <fc (/)and ^(J) -»/•„(/) so that 0 for any correlation iag. It means that the 
estimates are almost closed to the true values, which Is a sufficient performance. Note that 
the number of pilot symbols /^depends on the SNR. The higher the SNR is, the lower the 
number of N p is needed to reach a determineperformance. 

Estimator Architecture 

Fig. 2 represents a schematic architecture for an implementation of a RAKE receiver 
RECHV comprising a proposed channel estimator ESTIMOR using Implementing the 
algorithm (16). The estimator generates an estimate of p every slot. 

The RAKE receiver comprises: 

- a descrambler DESCRAMB, 

- a de-spreader CPICH_DESPREAD and DPCH_DESPREAD associated to each CPICH and 
DPCH channels respectively , in order to de-spread a received spreading signal, 

- the estimator ESTIMOR, and 

- a demultiplexer DEMUX In order to recover the two physical channels DPCCH and 
DPDCH. 

Note that a process of modulation of a transmitted signal is made with a pseudo- 
random sequence in order to limit the interferences between different base stations BS, this 
process, well known by the person skilled In the art, being called scrambling. At the receiver 
side, a de-scrambling needs to be achieved. 

The estimator ESTIMOR comprises: ' 
• Multipath coefficient estimate means CPICH_ESTTM that are adapted to calculate the 
estimate c^fa) of the multipath coefficient of step 2) from the CPICH channel's pilot 
symbol Scp^*(n), 

Factor estimate means BETA that are adapted to calculate the complex factor fi, 
A first Interpolator INTERPl on the CPICH branch, which Is enabled when the DPCH 
spreading factor Is equal to 512 i.e. when the CPICH symbol rate is lower than the DPCH 
one's, 
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- A second interpolator INTERP2 on the DPCCH branch, which is always enabled but not 
when the DPCH spreading factor is equalto either 256 or 512, i.e. when the DPCH 
symbol rate is lower than the CPICH one's. 

- Thus, when the spreading factor is equal to 256, none of the first or second 
interpolators are enabled, as it means that the CPICH and DPCH have the same symbol 
rate. 

In pracuce, in a preferred embodiment, nearest neighbor interpolation between two 
samples will be applied by the first INTERP1 and second INTERP2 interpolators, as it is 
shown to be a proper solution for this purpose. Indeed, it is the least complexity 
interpolation method while it yields a signal distortion due to such interpolation significantly 
below the noise level for practical ranges of Signal Interference (MAI) Noise Ratio SINR. 

In a non-limitative preferred embodiment, the estimator ESTIMOR further comprises: 

- Low-pass filtering means LPF1 and LPF2 for each channel CPICH and DPCCH 
respectively. 

The estimates of the cross- and auto-correlation (9) and (10) can be significantly 
improved by low-pass filtering successive estimates taken over several slot periods. Indeed, 
the beamforming weight vector nrhas a variation rate that depends on the variation rate of 
the pattYangles ft In practice path angles ft between the mobile terminal MT and the base 
stations BS, are likely to stay approximately constant over several slot periods i.e. there are 
only small angle variations. Slight beamforming adjustment may be needed although to 
track those small angle variations. Therefore p can be estimated and tracked by properly 
low-pass filtering successive estimates of the cross- and auto-correlation (9) and (10), e.g. 
with an exponential forgetting factor (which means that a smaller weighting Is applied on 
the older estimates), improving the noise suppression in the estimate. Note that as 
commonly known, any low-pass filter comprises a memory MEM to track the estimates. 

But, when the base station BS abruptly changes the beamforming weight vector w, 
the low-pass filters cant be used and usually the base station BS has to send a signaling 
MEM_RESET to the mobile terminal MT that provides a reset of the low-pass filter memory 
MEM. In case, such a signaling between the base station BS and a mobile terminal MT is not 
possible, the low-pass filters have to be removed. But whenever the low-pass filters are 
reset or removed, no signaling is needed to inform the mobile terminal MT toat transmit 
beamforming is used or not. Hence, said mobile terminal MT can keep estimating the factor 
/? independently of the presence of transmit beamforming, and when no transmit 
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beamforming is performed at the base station BS, fi will just be a real factor accounting for 
the transmit power offset between the CPICH and DPCH channels. 

In a non-limited preferred embodiment, the estimator ESTTMOR further comprises: 
5 - Two delay lines 2T 1 and Z -t(256/M) for the compensation of the delays Introduced by the 
low-pass filters LPF, the delay being knovyn and equal to L CPICH symbol periods. 

One can notice that the algorithm (16) required operations to generate an estimate 
of p every slot, which are 2N P products, 2{N P ^1) sums (N p being the number of pilot 

10 symbols), and 1 division In addition to the operations required to low-pass filter the cross- 
and auto-correlatlon estimates. In the case of the use of low-pass filters, the least 
complexity solution for this purpose consists of adopting, in a preferred embodiment, a first 
order IIR KALMAIM filter (with a suitable loop gain) as low pass filter, said filter being well- 
known by the person skill in the art The additional operations with this solution would be in 

15 general four products and two sums (i.e., two products and one sum per low-pass filter). 

In a non-limitative preferred embodiment, the estimator ESTIMOR further comprises: 

- A third interpolator INTERP3, which is enabled when the spreading factor of the DPCH 
channel is lower than 256, and which is adapted to interpolate one of the symbol rate of 

20 the estimate ft c^ ich (n) calculated at step 5) or of the symbol rate of the useful data of 

the DPDCH channel to match the other one. 

- A multiplier that multiplies the DPDCH mulOpath MP coefficient estimate to the data of 
said DPDCH channel, the output being sent to a combiner RECQV_COMB of the RAKE 
receiver. 

25 

Note that the structure in Fig. 2 is to be replicated for all RAKE receiver fingers. Each 
structure like the one in Fig. 2 has an output that is sent to the combiner RECEIV_COMB of 
the RAKE receiver in order to combine all the paths of a beam together. 

30 Conclusions 
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Before the optimization of a receiver, there is a phase of acquisition through a 
synchronization channel SCH, well known by the person of ordinary skill in the art, which 
consists in an estimate of the delays r p and the signal energy of the paths. 



According to the method of the Invention described above, we have Implicitly 
assumed that when the receiver can detect a path and estimate the associated delay on the 
CPICH through the synchronization channel SCH, a corresponding path at the same delay 
exists on the DPCH even though with a different complex amplitude. As we have seen, this 

5 assumption allowed the derivation of a DPCH channel estimation algorithm with low 

complexity, good performances and high flexibility, able to cope with the problem of channel 
estimation in the presence of transmit beamforming. 

It will be noted that under certain circumstances a path carrying significant signal 
energy on the CPICH can be detected at a certain delay by the synchronization channel SCH, 

10 while the corresponding DPCH path signal energy is too weak to be detected. Conversely, it 
may happen that a path that carries significant signal energy on the DPCH is actually 
undetectable on the CPICH, the acquisition phase being mute. The first case would lead to 
an estimate of p with a very small magnitude. Thus, the algorithm described above would 
still work properly and the RAKE receiver would simply drop the corresponding finger on the 

15 DPCH. The second case is much more critical because there would be no way to exploit the 
Information provided by the SCH and by the CPICH. Under these circumstances the 
estimator (13) cannot operate correctly and one should search for paths by directly using 
the second sub-channel DPCCH estimate, which is rather complex as it has been above- 
mentioned. 

20 

It Is to be understood that the present Invention Is not limited to the aforementioned 
embodiments and variations and modifications may be made without departing from the 
spirit and scope of the invention as defined in the appended daims. In the respect, the 
following dosing remarks are made. 
25 It is to be understood that the present invention is not limited to the aforementioned 

UMTS application. It can be use within any application of DS-CDMA using a receiver. 

It is to be understood that the method according to the present invention Is not 
limited to the aforementioned implementation! 

There are numerous ways of implementing functions of the method according to the 
30 invention by means of items of hardware or software, or both, provided that a single item of 
hardware or software can carries out several functions. It does not exdude that an assembly 
of items of hardware or software or both carry out a function. R>r example, the estimates of 
step 3) can be combined with the estimate of step 4), thus forming a single function without 
modifying the method of channel estimate with transmits beamforming in accordance with 
35 the invention. 
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Said hardware or software items can be implemented in several manners, such as by 
means of wired electronic circuits or by means of an integrated circuit that is suitable 
programmed respectively. The integrated circuit can be contained in a computer or in the 
receiver. In the second case, the receiver comprises means of the receiver adapted to make 
the de-scrambling, de-spreading, said means being hardware or software items as above 
stated, in the same way, the estimator within the receiver, comprises maximum likelihood 
estimate means adapted to make the ML estimates of step 1), first multipath coefficient 
estimate means adapted to make the multipath coefficient estimate of step 2), correlation 
estimate means adapted to make the correlations of step 3), complex factor estimate means 
that are adapted to make the factor complex estimate of step 4), and second multipath 
coefficient estimate means that are adapted to make the first sub-channel multipath 
coefficient estimate of step 5), as described previously, said means being hardware or 
software items as above stated. Of course, as mentioned above, said means can be an 
assembly of items of items of hardware or software or both carry out a function, or a single 
item carrying out several functions. 

The integrated circuit comprises a set of instructions. Thus, said set of instructions 
contained, for example, in a computer programming memory or in a receiver memory may 
cause the computer or the receiver to carry out the different steps of the estimate method. 

The set of instructions may be loaded into the programming memory by reading a 
data carrier such as, for example, a disk. A service provider can also make the set of 
instructions available via a communication network such as, for example, the Internet. 

Any reference sign In the following claims should not be construed as limiting the 
daim. It will be obvious that the use of the verb "to comprise- and its conjugations do not 
exclude the presence of any other steps or elements besides those defined in any claim. The 
word "a" or "an" preceding an element or step does not exclude the presence of a plurality 
of such elements or steps. 




1. A method for channel estimate with transmit beamforming, said channel comprising a 
first (CPICH) and second (DPCH) channels, said second channel (DPCH) comprising a 
first sub-channel (DPDCH) and a second sub-channel (DPCCH), characterized in that it 
comprises the steps of: 

• Model said first (CPICH) and second (DPCH) channels, said models being based on 
multipath (MP) components, a multipath component being characterized by a time- 
varying complex coefficient ( c „(0) and a delay (r p ), said time-varying complex 
coefficient having a first channel (CPICH) multipath coefficient (c wfcA (0 ), and a second 
channel (DPCH) multipath coefficient (c^fl) = P pCpepteh (t) ) depending on a complex 
factor (fi p ) due to the transmit beamforming, each p-th multipath component being 
processed as follows : 

1) Building a second channel (DPCH) and first channel (CPICH) instantaneous maximum 
likelihood (ML) channel multipath (MP) coefficients estimates {c dpdl (n), c cp[eb (n)); 

2) Building an estimate (c^w)) of the first channel (CPICH) multipath coefficient; 

3) Building estimates of a second channel (DPCH) and first channel (CPICH) cross- 
correlation (E{c dpell (l)^* (» -/)*}) and of a first channel (CPICH) auto-correlation 

(E^^OO^fa -/)*} ) from maximum likelihood (ML) estimates of step 2 at non- 
zero lag for noise suppression (x*0); 

4) Estimate said complex factor ( p p ) with said estimates In step 3; 

5) From steps 1 and 4 build a first physical channel (DPDCH) p-th multipath (MP) 
coefficient estimate {c 4pdeh (ri) = fi c^ri)). 

A method as claimed in dalm 1, characterized in that the second channel (DPCH) 
maximum likelihood (ML) channel multipath (MP) coefficient estimate (^(w) ) is only 

based on a maximum likelihood (ML) channel multipath (MP) coefficient estimate of the 
second sub-channel (DPCCH). 
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3. A method as claimed In claim 1, characterized in that the estimate of a second channel 
(DPCH) and first channel (CPICH) cross-correlation (E^Wc^ln -/)*} ) is an 
estimate of the second sub-channel (DPCCH) and first channel (CPICH) cross- 
correlation. 

5 

4. A method as claimed in dalms 1, characterized in that one of the instantaneous 
maximum likelihood (ML) channel multipath (MP) coefficients estimates of step 2 is 
interpolated to match a symbol rate of the other channel estimate. 

10 5. A method as claimed In claim 4, characterized in that an interpolation is done to the 
lower symbol rate between the first (CPICH) and second (DPCH) channels. 

6. A method as claimed In daim 1, characterized in that the second channel (DPCH) and 
first channel (CPICH) cross-correlation's estimate and the first channel (CPICH) auto- 

15 correlation's estimate are low-pass filtered. 

7. A method as daimed in claim 1, characterized in that a first channel (CPICH) multipath 
coefficient's estimate Is done by a linear prediction. 

20 8. A method as daimed in daim 1, characterized in that an estimate of said complex factor 
( p ) is achieved with a method of sample moments and ?«(*) ). 

9. A method as daimed in claim 1 or 8, characterized in that an estimate of said complex 

K 

factor ( p ) Is a weighed sum of estimates ( £ = A ) at lag I. 
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10. A method as daimed in claim 9, characterized In that an estimate of said complex factor 
( p ) Is limited to the first term of said weighed sum. 

11. A receiver utilizing said method as daimed in daim 1. 

12. A computer program product for a receiver, comprising a set of instructions, which, 
when loaded into said receiver, causes the receiver to carry out the method claimed in 
dalms 1 to 10. 
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13. A computer program product for a computer, comprising a set of instructions, which, 
when loaded into said computer, causes the computer to carry out the method claimed 
in claims 1 to 10. 

14. An estimator for channel estimate with transmit beamformlng, said channel comprising a 
first (CPICH) and second (DPCH) channels, said second channel (DPCH) comprising a 
first sub-channel (DPDCH) and a second sub-channel (DPCCH), characterized in that it 
comprises: 

- 1) Maximum likelihood estimate means that are adapted to calculate a second channel 
(DPCH) and first channel (CPICH) instantaneous maximum likelihood (ML) channel 
muitlpath (MP) coefficients estimates (c^(n), c^n)); 

• 2) First muitlpath coefficient estimate means (CPICH_ESTIM) that are adapted to 
calculate an estimate {c cplch (n)) of the first channel (CPICH) muitlpath coefficient; 

• 3) Correlations estimate means that are adapted to calculate a second channel (DPCH) 
and first channel (CPICH) cross-correlation (E{c^(0c^(» -/)*}) estimate and a 
first channel (CPICH) auto-correlation (E^Wc^/i-/)'}) estimate from said 
maximum likelihood (ML) estimates {c dpek (n) r c vleb (n)); 

■ 4) Complex factor estimate means that are adapted to calculate a complex factor ( p ) 
from said estimates of cross-correlation and auto-correlation; 

5) Second muitlpath coefficient estimate means that are adapted t calculate a first sub- 
channel (DPDCH) p*h muitlpath (MP) coefficient estimate (c^(w) = p c^iti) ) 

from the complex factor ( p ) and the estimate (c^ (») ) of the first channel (CPICH) 
multipath coefficient. 



15. A receiver comprising said estimator as claimed in daim 14. 
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Method for channel estimate with transmit beamforming 
ABSTRACT 

The present invention relates to a method for channel estimate with transmit 
beamforming, said channel comprising a first (CPICH) and second (DPCH) channels, said 
second channel (DPCH) comprising a first sub-channel (DPDCH) and a second sub-channel 
(DPCCH). It Is characterized in that it comprises the steps of: 

- Model said first (CPICH) and second (DPCH) channels, said models being based on 
multipath components, a multfpath component being characterized by a first channel 
(CPICH) multipath coefficient, and a second channel (DPCH) multipath coefficient 
depending on a complex factor (fi p ) due to the transmit beamforming, each multipath 
component being processed as follows : 

- 1) Building a second channel (DPCH) and first channel (CPICH) instantaneous maximum 
likelihood channel multipath coefficients estimates; 

- 2) Building an estimate (c„, fdi (/i)) of the first channel (CPICH) multipath coefficient; 

- 3) Building estimates of a second channel (DPCH) and first channel (CPICH) cross- 
correlation and of a first channel (CPICH) auto-correlation from step 2 at non-zero lag; 

- 4) Estimate said complex factor ( p ) from step 3; 

- 5) From steps 1 and 4 build a first physical channel (DPDCH) ^-th multipath coefficient 
estimate. 

Use: UMTS Receiver 

Reference: Fig. 2 
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